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ABSTRACT 

We present results from a blind, spectroscopic survey for z ^ 5.7 Lya -emitting galaxies using the 
Inamori Magellan Areal Camera and Spectrograph. A total of ~ 200 square arcminutes were observed in 
the COSMOS and LCIRS fields using a narrowband filter, which transmits between atmospheric emission 
lines at 8190A, and a mask with 100 longslits. This observing technique provides higher emission- 
line sensitivity than narrowband imaging and probes larger volumes than strong lensing. We find 170 
emission-line galaxies and identify their redshifts spectroscopically. We confirm three Lya emitting 
galaxies (LAEs), the first discovered using multislit-narrowband spectroscopy. Their line profiles are 
narrow, but fitted models suggest instrinsic, unattenuated widths ~ 400 km s~^ FWHM. The red wing 
of the line profiles present features consistent with galactic winds. The star formation rates of these 
galaxies are at least 5-7 M0 yr~^ and likely a factor of two higher. We estimate the number density 
of L > 5 X lO''^ ergs s'^ LAEs is 9.0^4^ x lO'^'h^Q Mpc'^ at redshift 5.7 and constrain the Schechter 
function parameters describing this population. Galaxies fainter than our detection limit may well be 
the primary source of ionizing photons at 2: ~ 6. We argue, however, that the break luminosity L*,li/q is 
not yet well constrained. If this break luminosity is near our detection limit, and somewhat lower than 
previous estimates, then the detected LAE population could be responsible for ionizing the intergalactic 
gas at redshift z ^ 6. We discuss the potential of multislit-narrowband spectroscopy for deeper emission- 
line surveys. 

Subject headings: galaxies: high-redshift - galaxies: luminosity function - techniques: spectroscopic — 
line: identification - line: profile 



1. INTRODUCTION 

Reionization of the intergalactic medium (IGM) marks 
a phase transition in the baryonic component of the uni- 
verse and quite likely heralds the emergence of the first 
galaxies that sustain star formation beyond a single burst 
(Oh & Haiman 2001). The appearance of complete Lya 
absorption troughs in quasar spectra at z> 6 shows that 
at least 1 in 10^ hydrogen atoms was still neutral then 
(Becker etal. 2001; Djorgovski etal. 2001; Fan etal. 2002; 
Songaila 2004; Fan et al. 2006), although much of the IGM 
appears to have been ionized hy z ^ 6.5 (Becker ct al. 2006, 
2007). The Thomson optical depth measured from the mi- 
crowave background polarization suggests an earlier start 
to reionization at Zr = III3 (Page et al. 2006; Spergel et al. 
2007). The Lya emission line, surprisingly, is proving very 
useful for identifying galaxies over this protracted period 
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of Reionization (lye et al. 2006; Stark etal. 2007b). Fre- 
quency shifting across the Stromgren spheres of the sources 
allows a fraction of the emission line photons to propagate 
through an IGM with a considerable neutral fraction. Un- 
derstanding the intrinsic evolution of Lya emitters (LAEs) 
builds a foundation for using the apparent evolution of the 
LAE luminosity function during reionization to measure 
the changing neutral fraction (Malhotra & Rhoads 2004; 
Kashikawa et al. 2006). 

High-redshift LAEs may provide the first glimpse of 
primeval galaxies, owing to the high Lya equivalent- width 
of young, metal-poor galaxies (Malhotra & Rhoads 2002; 
Dawson etal. 2004; Schaerer 2007). The stellar mass den- 
sity at z ~ 5 (Stark etal. 2007a), the presence of carbon 
in the IGM at z = 5.7 (Simcoe 2006; Ryan- Weber, Pet- 
tini, & Madau 2006), and the steep faint-end slope of the 
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i-dropouts (Bouwens et al. 2007) can all be interpreted as 
evidence for a significant population of as yet undetected 
dwarf galaxies at high redshift. The properties of these 
low-mass galaxies during and just after reionization are 
cosmologically important. Finding them and measuring 
their properties will continue to be of interest well into 
the next decade. 

Finding high-redshift galaxies remains largely a contrast 
problem with the sky. Selection by Lya emission provides 
higher sensitivity than continuum break identification for 
populations with large line equivalent widths, like young, 
metal-poor galaxies. Over the past decade, several tens 
of LAEs have been confirmed near redshift 6 among ob- 
jects found via narrowband imaging (Hu et al. 2002; Cuby 
et al. 2003; Ajiki et al. 2003; Kodaira et al. 2003; Hu et al. 
2004; Rhoads et al. 2004; Ajiki et al. 2004; Taniguchi et al. 
2005; Shimasaku etal. 2006; Kashikawa etal. 2006; Mu- 
rayama etal. 2007; Westra etal. 2005, 2006; Ouchi etal. 
2007). The redshifts of nearly all of these galaxies shift the 
Lya line into one of the bandgaps in the OH (hydroxyl) 
spectrum of the terrestrial sky. The dark regions at 8200A 
and 9200A are the widest 150 A) and are useful for 
studying the end of reionization. Lyman-alpha emission 
at z = 5.7 and 6.5 redshifts into these bands. 

Spectroscopic surveys also benefit from these low back- 
ground regions and offer higher sensitivity than narrow- 
band imaging. Dispersion of the background improves 
line-to-sky contrast, although a smaller solid angle of sky is 
covered. One prolific spectroscopic approach has been tar- 
getting strongly lensed LAEs. Probing comparatively tiny 
volumes, three very faint LAEs have been confirmed with 
this technique at z 5 — 6 (Ellis etal. 2001; Santos etal. 
2004). Using the Hubble Space Telescope to get above 
the atmospheric OH emission has also produced 3 LAEs 
at z > 5 (Pirzkal etal. 2007). Serendipitous discoveries 
have also been made in the "sky region" of longslit spec- 
tra (Sawicki etal. 2007; Stern etal. 2007). In this paper, 
we present the first detections of LAEs with a technique 
- multislit-narrowband spectroscopy (MNS) - designed to 
optimize the sensitivity of ground-based telescopes to high- 
redshift LAEs. Previous surveys with MNS failed to find 
any LAEs (e.g. Crampton & Lilly 1999; Martin & Sawicki 
2004; Tran et al. 2004). 

Observations were made with the Inamori Magellan 
Areal Camera and Spectrograph (IMACS, Dressier etal. 
2006) on the Baade Telescope, which offered a wide field 
of view. Longslit spectra were obtained in the dark re- 
gion at 8200A. A custom blocking filter and "Venetian 
blinds" mask multiplexed the effective area one-hundred 
fold. We carried out four blind searches: two in the COS- 
MOS (Cosmic Evolution Survey) field and two in a LCIRS 
(Las Campanas Infrared Survey) field. The total survey 
area of 200 square arcminutes is significantly larger than 
that attempted previously with similar techniques. 

We describe the survey in §2 and explain how we found 
emission-line objects and identified their redshifts in §3. 
Some readers may wish to skip directly to §4, where we 
discuss the star formation and feedback properties of the 
LAEs, estimate their space density, and explore the con- 
tribution of the LAE population to the ionization of the 
IGM at z = 5.7. We summarize the advantages and disad- 
vantages of MNS in §5, including a direct comparison to 



planned instrumentation for the next generation of tele- 
scopes. Section [6] summarizes our conclusions. Proper- 
ties of the foreground line-emitters will be discussed in a 
separate paper. We adopt a cosmology with fl„i = 0.3, 
flA = 0.7, Ho = 70 km s'^ Mpc'^, and flbh^ = 0.02 
throughout this paper. The angular scale at redshift 5.7 
is then 5.87 kpc/"; the universe is 0.98 Gyr old; and the 
time between redshift 6.5 and 5.7 is 151 Myr. 

2. OBSERVATIONS AND DATA ANALYSIS 

Four types of data were used in our emission-line sur- 
vey. We began with the blind, emission-line survey with 
the IMACS instrument on the Baade telescope. Broad 
bandpass spectroscpy of emission-line objects was then ob- 
tained to identify the redshift of the line via the presence 
(or absence) of other lines. A subsample of LAE candi- 
dates was then followed up with higher-resolution spec- 
troscopy to determine the shape of the line profile. Broad- 
band imaging of our two fields by the COSMOS (Scoville 
et al. 2007; Taniguchi et al. 2007) and the LCIRS (Marzke 
etal. 1999; Chen etal. 2002) supplement our new spec- 
troscopic observations. The photometric redshifts from 
COSMOS (Mobasher et al. 2007) were particularly use- 
ful because the association of a line-selected object with a 
continuum source is often ambiguous. 

2.1. Multislit, Narrowband Spectroscopic Search 

Multislit narrowband spectroscopy remains a relatively 
new technique, and we illustrate the filter transmission 
in Figures [1] and the mask design in Figure [5) Our cus- 
tom blocking filter fabricated by Barr Associates passes 
a 121 A FWHM bandpass at 8188 A. The transmission 
is over 90% across this band and drops sharply to less 
than 5% transmission outside the A8104 - A8273 band- 
pass. We placed this filter in the parallel beam directly in 
front of the 200 lines mm~^ grism and recorded data with 
the short, f/2, camera. With a dispersion of 2.0 A pix~^, 
we fit spectra from 100 longslits (numbered 9 to 108) on 
the eight 2k by 4k thinned, SITe detectors. At low disper- 
sion, the angle between the first-order spectrum and the 
zero-order image of the mask is smaller than the detector, 
so a zero-order image of the mask also falls on the de- 
tector. This zero-order image is easily recognized via the 
pattern of the multislit mask, and this region is masked 
and ignored in the analysis. The 1'.'5 wide slits subtend a 

geometric area of 55.3 □ per pointing within the 27f50 di- 
ameter field of the camera. With this setup, emission lines 
from diffuse sources have a width of 12.8 ± 0.3 A FWHM. 
The spectral resolution is typically a factor of two bet- 
ter for spatially-unresolved objects, as determined by the 
atmospheric seeing. 

Table [1] summarizes the search observations. Two masks 
were exposed in each of two fields. The masks were offset 
by half the slit separation. Exposure times ranged from 
6 to 11 hours, and the image quality was generally 0'.'8 
FWHM or better. Image quality was good across the de- 
tector during the 2005 March and May runs, which used 
the field corrector and atmospheric dispersion compen- 
sator. The 2004 April run took place prior to installa- 
tion of the ADC/corrector, and the image quality was fur- 
ther compromised by an un-even focus across the detector 
caused by a misaligned tertiary mirror. 
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We used techniques standard to optical imaging to re- 
move the instrumental signature from the data. The band- 
pass is narrow enough to ignore the wavelength depen- 
dence of the detector response; and we corrected for pixel- 
to-pixel response variations with direct images of the flat- 
field screen taken through the NB8190 filter (in the filter 
wheel). Stacked images were constructed using cosmic ray 
masks. Faint, sky lines are present even in the relatively 
clean, dark bandpass at 8200A. We subtracted this sky 
background after computing it in a local aperture along the 
slit direction. Subtraction of this sky model revealed the 
objects. These procedures were implemented with a set 
of custom scripts, implemented with PerlDlQ and IRAF 
routines 

We used the location of each object along the slit and 
the transformation from mask to sky coordinates, kindly 
provided by K. Clardy, to determine object positions in the 
spatial direction. We tested the accuracy of the relative as- 
trometry using bright, foreground objects that fell on slits 
and the broadband images described below. The position 
errors along the slit were normally distributed with stan- 
dard deviation of 0'.'24 for the 2004 data, slightly less for 
the 2005 observations. The same bright objects present 
offsets perpendicular to the slit of ± 0'.'36 (Itr result). We 
emphasize, however, the degeneracy between the position 
of an object in the slit and the wavelength of the emission 
line that we measure. Half a slit width is ± 0'.'75, and 
offsets approaching this limit lead to wavelength errors 
< 7. 5 A. The accuracy of the zero-point of the astrometric 
transformation was estimated by comparing the coordi- 
nates of the mask alignment stars, taken from the USNO 
catalog, to their coordinates in the Subaru images of the 
COSMOS field. After rejection of two stars, we find a sys- 
tematic offset of 0'.'40 W and and 0'.'29 S of the Subaru 
positions using four stars. 

We inspected the stacked, sky-subtracted frames from 
each of the four narrowband, multislit observations for 
emission lines. Line detection was then automated with 
SExtractor (version 2.3.2, Bertin & Arnouts f996), requir- 
ing detection over at least 9-contiguous pixels, correspond- 
ing to a spatial width > 0'.'6 and a line width > 6.4 A. 
The SExtractor software recovered our initial, by-eye list of 
sources; additional candidates were inspected and found to 
be associated with incomplete masking of chip edges, zero- 
order images, and cosmic-rays. The SExtractor software 
worked well for finding the emission-line spectra because 
the signal from a line is very similar to that from an object 
in a direct image. The photometry from SExtractor is not 
appropriate for our survey however, because sky and con- 
tinuum subtraction needed to be handled with techniques 
specific to spectroscopic data. 

We measured the flux of each emission line in a fixed 
apertures using IRAF. This line photometry was per- 
formed on the images prior to sky subtraction using aper- 
tures defined on the sky-subtracted frames. Because faint. 



sky lines are present even in this relatively clean window, 
care had to be taken to measure the average sky counts 
from exactly the same bandpass as the line aperture. Con- 
tinuum was measured in the same spatial aperture as the 
line but was usually undetected. A flux calibration was de- 
rived from observations of spectrophotometric standards 
taken during each run. The throughput during the April 
2004 run was 28% better than that during the spring 2005 
runs due to a recent re-aluminization of the telescope mir- 
ror. The uncertainty in this calibration is 5-10% (except 
for the 2005 May search), and our photometric errors are 
normally dominated by the random pixel-to-pixel varia- 
tions in the background under the object and/or the Pois- 
son noise from the total counts. The May run was gen- 
erally clear but not photometric, and we adopt the flux 
calibration derived 2 months earlier with the identical in- 
strumental setup. 

2.2. Broad Bandpass Spectroscopy 

We obtained spectra of emission-line candidates with 
IMACS, a low resolution grism (usually 200 mm~^), the 
GG455 blocking filter, and multislit masks with 1'.'5 slits. 
Masks for the lOH field were observed in 2005 May (4.5 
hrs) and 2007 March (5.0 hrs). Masks in the 15H field were 
observed in 2004 July (7.0 hrs), 2005 June (2.0 hrs), and 
2006 June (3.08 hrs with the 150 mm"-'^ grism). The wave- 
length coverage for a typical source was 5500 A to 9500 A. 
The C0SM0q3 data reduction package was used to stack 
frames, remove cosmic rays, subtract sky, calibrate wave- 
length, and extract spectra. The 2007 March and 2004 
July observations of the lOH and 15H fields, respectively, 
are significantly more sensitive than the other follow-up 
observations. 

To date, low-resolution spectra have been attempted for 
95% of the emission-line objects in both the lOH and 15H 
fields with a recovery rate of 80% and 95%, respectively. 
Relative astrometric errors were of particular concern for 
pre-corrector discoveries in the 2004 data, but we found no 
tendency for unrecovered objects to be on the periphery 
of the field. 



We 



2.3. Higher Resolution Spectroscopy 
obtained higher resolution spectra of 



Lya candidates in the lOH Field with the IMACS long, 
f/4 camera and the 1200 mm^^ grating in 2006 Febru- 
ary. Line emission from sources filling the 1'.'5 slits has a 
FWHM of 94 km s~^ ; and the resolution for unresolved 
objects is 50 to 63 km s^^ . The two masks included only 
a subsample of the candidates passing the f/2 follow-up. 
Exposure times of 8.24 hrs and 8.75 hrs were sufficient un- 
der clear conditions to detect most targets. Some data for 
the 15H Field were obtained with the same configuration 
in 2006 February (6.0 hrs and some cloud cover) and 2006 
June (3.5 hrs with cirrus; and 5.0 hrs under clear skies). 



^ The Perl Data Language (PDL) has been developed by K. Glazebrook, J. Brinchmann, J. Cerney, C. DeForest, D. Hunt, T. Jenness, T. Luka, 
R. Schwebel, and C. Soellcr and provides a high-level numerical functionality in the perl scripting language (Glazebrook & Economou, 1997). 
It can be obtained from http:\\ pdl.perl.org. 

^ The Image Reduction and Analysis Facility, a general purpose software system for the reduction and analysis of astronomical data. IRAF 
is written and supported by the National Optical Astronomy Observatories, which is operated by the Association of Universities for Research 
in Astronomy, Inc. under cooperative agreement with the National Science Foundation. 

* The Carnegie Observatories System for MultiObject Spectroscopy was created by A. Oemler, K. Clardy, D. Kelson, and G. Walth; it is 
available at http://www.ociw.edu/Code/cosmos . 
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Resolved spectra in the 15H field were obtained at Keck 
in 2006 May using the 831 mm~^ grating on LRIS (Oke 
etal. 1995). Under thin clouds, integration times of 1.5 to 
3 hrs yielded detections of most sources. Through the 1'.'5 
slits, telluric emission-lines had a FWHM of 143 km s^^ , 
and spectra of unresolved sources had spectral resolution 
of 76 to 95 km s'^ . 

2.4. Continuum Imaging 

Deep, broadband imaging complements our emission- 
line survey because it is more sensitive to continuum emis- 
sion. Photometry can confirm redshifts via identification 
of the Lyman break, measure star formation rates via the 
strength of the rest-frame UV continuum, and constrain 
age via the Balmer break. Scattering of Lya photons by 
intergalactic hydrogen at 4.0 < z < 5.7 attenuates the 
continuum between the Lyman break and the Lya line, 
allowing only ^ 20% transmission (Fan etal. 2002). At 
redshift 5.7, the Lya line falls in the i (and i or /) band, 
so imaging in z (or IR) filters provides UV (optical) con- 
tinuum measurements. Much of the r (and r and R) 
bandpass lies below the Lyman limit. Only filters pass- 
ing light shortward of the Lyman limit provide good veto 
imaging for high-redshift candidates. 

2.4.1. 1 OH Field 

Taniguchi etal. 2007 present Subaru/Suprime-Cam 
imaging of the COSMOS field in the UBgVrIz broad- 
bands and one narrowband NB816. The imaging redward 
of our emission lines reaches z = 25.9 (Taniguchi et al. 
2007; 2"diameter apertures; 3a), which is deeper than the 
z sa 25.0 (S. Gwyn website)0 imaging from the Canada 
France Hawaii Telescope Legacy Survey (CFHTLS-Field 
D2). Most of the emission-line objects that we find are 
not bright enough, > 5 x 10^^^ erg s^^ cm^^, to be de- 
tected in the COSMOS z' band. 

Assuming flat or redder SEDs (spectral energy distribu- 
tions) in /y, objects not detected in the Subaru z imaging 
will not be detected in either the Subaru i imaging or the 
HST ACS F814W observations. The latter are slightly 
more sensitive than the Subaru imaging for sources less 
than I'.'O in diameter. The bandpass of our 8190 A fil- 
ter lies just within the red limit of these filters, so these 
broadband filters sample mainly attenuated continuum. 
The Bg'v bands from COSMOS/Suprime-Cam are below 
the Lyman limit of z = 5.7 objects, so a detection would 
veto identification of the object as a high-redshift galaxy. 

2.4.2. 15H Field 

Far less imaging is available for the 15H Field. The 
z imaging from LCIRS is significantly deeper than the 
infrared imaging. It reaches z — 23.5, equivalent to 
z'ab = 24.0. 

We obtained a new, deep image in a 3000 A wide filter, 
WB48-78, in 2005 March. A total of 3.5 hr of integration 
was obtained with the IMACS short camera on Magellan. 
The blue half of this bandpass lies longward of the Lyman 
limit, but the red half includes the Lya and Ly/3 atten- 
uated continuum. Null detections in the WB48-78 image 



can therefore confirm a continuum break, but a detection 
does not necessarily rule out a continuum break. 

3. RESULTS 

3.1. Detection of Emission Lines 

Over cosmological distances, the HI opacity of the IGM 
blankets the continuum from the Lyman limit to the 
Lya line, effectively shifting the continuum break up to 
the Lya line (Madau 1995). Any continuum blueward 
of a redshift 5.7 Lya line must be weak due to Lya 
scattering along the line-of-sight. None of the line-plus- 
continuum spectra present continuum breaks across the 
emission line, so we excluded all the line-plus-continuum 
sources from the preliminary list of 136 Lya candidates. 
Few of the line-only objects in the lOH field were detected 
in the Subaru narrowband image confirming our argument 
that spectroscopic surveys offer unique emission-line sen- 
sitivity. Table [5] shows the object yield for each of the four 
search runs, distinguishing line emitters with and without 
continuum in their discovery spectra. 

3.2. Confirmed LAEs 

We have confirmed three LAEs among the emission- 
fine galaxies: MSDM80-f3 in the COSMOS field and 
MSDM29.5-h5 and MSDM71-5 in the 15H field. Confirma- 
tion requires: (1) the broad bandpass spectrum presents 
no spectral lines between ~ 5500 and 9500A other than 
the discovery line near 8200A and (2) the profile of this 
line, at higher dispersion, is skewed to red wavelengths, 
consistent with transmission through high-redshift IGM. 
In this section, we present spectra of these three objects, 
put limits on the number of LAE candidates, and derive 
the survey completeness. Properties of the LAEs are fur- 
ther discussed in § 14.11 

Lyman-alpha-emitter MSDM29.5+5 from the 15H field 
shows the classic red-skewed line asymmetry of a LAE in 
the high-resolution, Keck spectrum; see panel (c) of Fig- 
ure H The line width is 235 km s'^ FWHM. The low 
resolution spectrum, panels (a) and (b), is not sufficient 
to recognize the asymmetry. The line in the 200-1 spec- 
trum is narrower than the overlaid sky line, indicating a 
source size considerably smaller than the slit. The atten- 
uated Lya flux is ~ 1.8 ± 0.2 x 10^^*" erg s^^ cm^^ in 
the 2005 May discovery spectrum. The systematic errors 
in the flux calibration may be more uncertain than this 
statistical error due to the observing conditions described 
in § m The observed-frame equivalent width is at least 
11 A based on the 3cr upper limit on the continuum in the 
search spectrum. This object is outside the fleld of the 
LCIRS image. 

Figure g] shows LAE MSDM80+3, discovered in the 
2004 search of the lOH field. The broad bandpass spec- 
trum shows a slight red asymmetry of questionable signifi- 
cance. The higher resolution spectrum shows multiple line 
components spread over about 400 km s^^ FWHM. The 
Lya flux is 1.4±0.2 x 10^^'' erg s^^ cm^^. The 3cr limiting 
z flux (for 2" diameter apertures) gives a lower limit on 
the observed-frame equivalent width of 144 A. 

We believe MSDM80+3 is the same object as 
#55 detected by Muramyama etal. (2007) using 



^ |http://www3.Ccidc-ccda.hia-iha.nrc-cnrc.gc.ca/community/CFHTLS-SG/docs/cflitls.html| 
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Subaru/Suprime-Cam narrowband imaging. The dechna- 
tions agree to 0'.'2. Their right ascension is nearly I'.'O 
west of ours, consistent with the zeropoint shift we found 
relative to the Subaru image in § [2l They estimate a sig- 
nificantly higher line flux of 2.8±0.4x 10~^^ erg s^^ cm^^. 
The discrepancy could result from attenuation of the 
source by our slit. Line fluxes from narrowband imag- 
ing are, however, subject to uncertainties about redshift, 
filter transmission, and continuum level. Since no contin- 
uum is detected from MSDM80-I-3 in broadband i or z 
images (Taniguchi et al. 2007) and we find a wavelength 
near the center of the filter bandpass, we see no obvious 
reason for the imaging result to be in error. We quantify 
our slit losses on a statistical basis in § 13.41 and include 
them in modeling the number counts in § 14.21 

Another line-emitter, MSDM71-5, marginally meets our 
criteria to be labeled "confirmed." The shape of the line- 
profile leaves some ambiguity. Figure [5] shows follow-up 
spectra of MSDM7f-5 taken with the 200-1 and the 150-1 
grisms. Even at such low resolution, both of these spectra 
hint at a broad, red wing. Curiously, our high resolution 
Keck, LRIS spectrum shows a narrow line, ft is highly 
unlikely that different objects were observed with the two 
telescopes because the line appears at the same wavelength 
in all 3 spectra. The B and z LCIRS images show no ob- 
jects near the position of MSDM71-5. We detect an object 
slightly less than 1'.'5 to the southeast in our deeper WB48- 
78 image; but, as explained in § 12.4. 2[ a weak detection in 
this band does not veto a line identification as Lya . 

The coordinates, approximate wavelength, and line flux 
of these confirmed LAEs are given in Table [H We also 
tabulate the properties of emission line objects for which 
identification of the line as Lya has not been ruled out, 
as we explain in the next section. 

3.3. Additional LAEs? 

We now ask whether our selection criteria may have 
missed some LAEs. Our follow-up observations were un- 
usually complete. Low resolution, follow-up spectroscopy 
has been attempted for all but three line-only objects 
in each field. Should we have rejected all the line-plus- 
continuum objects? Many line-plus-continuum objects 
were in fact used to fill our follow-up masks. All their 
spectra obtained to date have been identified as foreground 
objects (see Table [5]), so our initial continuum cut appears 
to be valid. The only sources that allow further specu- 
lation are the un-recovered targets and the galaxies that 
present a single emission line, i.e. just the discovery line 
near 8200A, in their broad bandpass follow-up spectrum. 

Table [2] reveals that in the lOH and 15H fields, 13 and 7 
candidates, respectively, were not recovered by follow-up 
observations. Some candidates with line fluxes near the 
detection limit, derived in § 13. 4| will inevitably turn out 
to be spurious. Review of the stacked, search data, led 
us to believe the majority of such sources were genuine, 
however. We attribute some failed recoveries, particularly 
for the 2004 sources, to the relative astrometric errors de- 
scribed in ij 12.21 The fraction of LAEs among these objects 
is presumably as small as it is for the full sample, so we 
estimate 0-1 LAEs were missed. Our follow-up observa- 
tions are also biased against detecting the faintest LAEs 
in our sample. Time and weather limitations resulted in 



broad bandpass spectroscopy with half the exposure time 
of the initial search. These objects may include a higher, 
but unknown, fraction of LAEs. 

Another set of 14 candidates (7 in the lOH field and 
7 in the 15H field) are listed in Table [3l These objects 
present a single emission-line, at ~ 8200 A, in their broad 
bandpass spectrum. Higher resolution spectra, had they 
been obtained for these objects, would differentiate be- 
tween LAEs, foreground [Oil] emitters at z = 1.2, and 
extremely faint [OIII] 5007 emitters - i.e. those with un- 
detectable [OIII] 4959 and H/3 lines. Among the objects 
observed at higher-dispersion, the yield of LAEs to [Oil] 
emitters varied widely among fields, making it difficult to 
predict the LAE fraction among the single- line sources. 
We gained additional insight from photometric rcdshifts 
from COSMOS (Mobasher et al. 2007) in our lOH field. 
The continuum emission from our candidates is typically 
fainter than the i < 25 completeness limit. We did investi- 
gate the redshifts of foreground galaxies near the sightline 
to each line emitter, however, with the idea that associ- 
ated satellite galaxies or HII regions might be the source 
of the detected line emission. We find foreground galax- 
ies surprisingly close to many of the single-emission-line 
galaxies, and their photometric redshifts often allow [Oil] 
emission to be shifted into our search bandpass. For read- 
ers interested in the exact accounting, sections 13.3.11 and 
I3.3.2l review the single-emission-line objects individually in 
the lOH and 15H fields, respectively. 

3.3.1. Single-Emission-Line Candidates in the 1 OH Field 

The low-resolution line profiles of the seven single-line 
candidates in the lOH Field are shown in Figure [HI The 
line profiles of MSDM52-H8, MSDM17.5-h8, and MSDM13- 
3 present significant skew. The spacing of the two peaks 
in the MSDM52-f 8 spectrum is consistent with [Oil] emis- 
sion; however, the bluer of the two emission peaks is 
stronger, which would imply a high electron density if the 
line is indeed [Oil]. Electron densities in typical HII re- 
gions are lower than this, so this line can only be [Oil] if 
it is an AGN (Active Galactic Nucleus), which will rarely 
be the case. In contrast, the profile shape of MSDM13-3 is 
consistent with a blended [OH] doublet in the low-density 
limit. These spectra lack the resolution required to reveal 
the asymmetry of many LAEs, such as MSDM29. 5-1-5 in 
Figure [21 so a larger fraction of the single- line sources may 
actually have asymmetric line profiles. 

To identify more interlopers, we searched the photomet- 
ric redshift catalog of Mobasher et al. (2007) near the po- 
sitions of our emission-line objects. The density of objects 
in the Subaru (UBgVrIz) imaging is such that a foreground 
galaxy brighter than the i=25 detection limit lies within 
4'.'5 (or 9'.'0) of a randomly chosen position 51% (or 92%) 
of the time. The projected distance to the nearest red- 
shift 1.2, [OII]-emitting galaxy is larger. In the Mobasher 
etal. catalog, we find a galaxy consistent with being at 
redshift 1.2 within 4'.'5 (or 9'.'0) of a random sky position 
only 3.5% (or 11%) of the time. In light of these statistics, 
if the single-line objects are primarily LAEs, then very few 
(if any) will be projected near a galaxy with a photometric 
redshift near 1.2. 

As expected, no foreground continuum sources are de- 
tected near our one confirmed LAE in the lOH field, 
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MSDM80+3. The closest galaxy whose photometric red- 
shift is consistent with [Oil] emission in our filter lies 
9'.'9 to the northeast. Our position for MSDM64-5 is 2'.'1 
northwest of the center of a foreground spiral galaxy with 
Zphot = 0.90. The 95% confidence interval, 0.82 < z < 
1.02, allows H7 (and plausibly a faint [OIII] 4363 line) to 
fall in our search band. The problem with this associa- 
tion is that H/3 and [Oil] were not detected in the broad 
bandpass spectrum. 

The other six single-line galaxies in the lOH field are 
likely foreground galaxies. All of them lie less than 3'.'5 
from a galaxy plausibly at redshift ~ 1.2. The associ- 
ations are particularly close, in both angular separation 
and probable redshift, for MSDM30.5-8, MSDM17.5-t-8, 
and MSDM17-1. The separation is larger, 2'.'1 (17.1 kpc 
at z = 1.2), for MSDM52+8, but much of this discrep- 
ancy is in the direction of our astromctric (zero-point) er- 
ror. We argued that MSDM13-3 is likely an [Oil] emitter 
based on its line profile. Both MSDM13-3 and MSDM50-7 
are very close to faint, i > 25, galaxies, which are known 
to be foreground owing to B-band detections. Theses off- 
sets are small enough to be astrometric errors, and we 
tentatively assign our emission-line detection to the faint 
objects. These i > 25 galaxies may be associated with the 
brighter, cataloged galaxies (near redshift 1.2) a few arc- 
seconds away. Further work will be required to understand 
their physical nature. 

3.3.2. Single- Emission- Line Candidats in the 15H Field 

In addition to MSDM29.5-h5 and MSDM71-5, we iden- 
tified 7 single emission-line sources in the 15H field. Their 
f/2 spectra are shown in Figure [71 The line profiles of 
MSDM57+4, and MSDM70-M show some sign of asym- 
metry at low resolution. We interpret the blue skew of 
the MSDM70-hl profile as evidence of the [Oil] doublet, 
an identification consistent with the source being the ob- 
ject in the (LCIRS) B band image, located less than 1'.'5 
from our estimated position. The only other single-line 
object with a B-band candidate this close to our position 
is MSDM89-2, which we also exclude from the candidate 
LAEs on the basis of a probable B-band detection. 

All the single-line candidates show at least a faint 
smudge in the WB48-78 image within 1'.'5 of our position. 
If these faint objects are the counterparts of our LAEs, 
they should be brighter in the z band. These galaxies 
are not detected in the LCIRS-z images, but the depth 
is not sufficient to rule out a continuum break across the 
line. It is only the high frequency of plausible [Oil] inter- 
lopers in the lOH field that strongly suggests many will 
turn out to be foreground interlopers. The extreme depth 
of the WB48-78 image does demonstrate that plausible 
continuum counterparts can be found at offsets from the 
emission-line objects that are significantly smaller than 
those from bright objects with photometric-redshifts in 
the lOH field. This result strengthens our suspicion that 
the objects near redshift 1.2 that are several arcseconds 
from our emission-line position are not the same source, 
although the two objects may be associated (i.e. at the 
same distance). 

3.4. Completeness and Effective Area 

To determine completeness levels, we added objects of 
various brightness to the search frames and then checked 



whether SExtractor recovered them. The intensity of an 
object was modeled by a gaussian profile ranging in size 
from 0'.'45 to 0'.'90 spatially and linewidth 6A to 12A. We 
tuned SExtractor to find most of our by-eye sources but 
not to go too much fainter, which yields many spurious 
sources. This tuning requires that an object have at least 
10 pixels above 1.5 times the background rms. SExtractor 
found 90% of our by-eye sources and most of the 'missing' 
10% can be explained as being near image artifacts that 
confuse SEextractor. This tuning of SExtractor allowed 
us to recover the simulated sources in a way that was both 
automated and very close to what the human eye finds. 

Results of the Monte Carlo simulations are shown in Fig- 
ure [8] Most of our sources were unresolved spatially, so we 
interpolate to the results for 0'.'7 seeing. The recovery rate 
of our survey is 50% at a line flux of 6 x 10~^^ erg s~^ cm~^ 
for the most sensitive search, the 2004 observations of the 
15II field, degrading to 9 x 10~^® erg s^^ cm^^ for the 
least sensitive search, the 2005 15H field. 

Attenuation by the slit and masked areas such as chip 
edges, the zero-order image of the mask, and bright fore- 
ground objects affect the effective area of our survey. The 
Monte Carlo simulations quantify how the slit length that 
we can search is reduced by masked objects. For a fixed slit 
width of 1'.'5, this effective slit length per IMACS pointing 

ranged from 45.1340 in the 15II 2005 search to 45.589 

□ in the lOH 2005 search. The effective area over which 
a particular source could have been detected depends on 
its brightness. Emission-line objects near our detection 
threshold will only be detected if a slit is centered on them, 
while extremely bright sources could be highly attenuated 
and still detected. These slit losses were modeled using a 
Gaussian brightness profile as illustrated in Figure [9l 

We modeled the effective slit length and completeness 
separately for each of our four search runs. The calcula- 
tions in the next section use the results. 

4. THE LAE GALAXY POPULATION 

In this section, we discuss the properties of the LAEs, 
compute their space density, and find the families of 
Schechter function parameters that can describe this 
galaxy population. These descriptions of the LAE pop- 
ulation are used to estimate the production rate of Lyman 
continuum photons, which is compared to the ionization 
requirement of the IGM at redshift 5.7. 

4.1. Properties of LAEs 

The Lya fluxes that we measured correspond to lu- 
minosities of logLya (erg s'^) = 42.70, 42.81, and 42.89 
for MSDM80-f3, MSDM29.5-h5, and MSDM71-5, respec- 
tively. These LAEs are fainter than the luminosity break 
(i.e. L*) fitted to a larger sample of LAEs in the Sub- 
aru Deep Field (Shimasaku etal. 2006). Other surveys 
that have confirmed LAEs this faint at z ~ 5.7 include: 
GRAPES - the Grism ACS Program for Extragalactic Sci- 
ence (Pirzkal et al. 2007) in the Hubble Ultra Deep Field 
(HUDF, Beckwith et al. 2006), the Gemini Deep Deep Sur- 
vey (GDDS) ~ which obtained 30 hr spectrosopic observa- 
tions of i-dropouts in the HUDF (Stanway etal. 2006), 
a survey of cluster caustics for lensed LAEs (Ellis et al. 
2001), and the narrowband imaging surveys of Hu etal. 
(2004) and Shimasaku et al. (2006) with Subaru. 
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We estiraate star formation rates for the Magellan LAEs 
from their line luminosities. Carrying the amount of Lya 
attenuation, fiya, along as a parameter, the Kennicutt 
(1998) conversion from Ha luminosity to SFR and the 
Lya to Ha flux ratio from Case B recombination (Brock- 
lehurst 1971) yield star-formation rate estimates 

SFR = 1.0/^;„ Mq yr-i(Li^„/1.05 x 10^^ g^gs s-i)(l) 
This expression assumes continuous star formation, so- 
lar metallicity, and a stellar initial mass function with 
dlogN/dlogM = -2.35 from 0.1 to 100 Mq. The es- 
timated SFRs of the LAEs are 4.8/;~^, 6.1/;^^, and 
7Af~l Mq yr"^ Variations in the ratio L{Lya) / L^iJJV) 
are large among LAEs, but star formation rates derived 
from the UV continuum are typically a factor of a few 
times higher than those estimated from Lya (e.g. Pirzkal 
etal. 2007). The faint LAE population is very blue 
(Stanway etal. 2006; Pirzkal etal. 2007), so the upward 
corrections to these UV luminosities for extinction may 
be small. Our guess for the new LAEs is an attenu- 
ation fiya ^ 0.25 — 0.5, which corresponds to SFRs of 
-10-30 Mq yr-i. 

Broadband photometry of the COSMOS field (and the 
LCIRS field) is not currently deep enough to detect the 
LAEs from the Magellan MNS survey. Rest-frame UV-to- 
optical colors are available for nine LAEs in the HUDF at 
slightly lower redshift 4.0 < z < 5.7. Since the Lya prop- 
erties of the two samples are similar, the stellar masses 
estimated by Pirzkal et al. (2007) for the former, 10^ to 
10^ Mq, may be indicative of the masses of the Magellan 
MNS LAEs. The LAEs are not resolved in good seeing, so 
the size of the emitting region is less than about 5 kpc. 

The shapes of the Lya line profiles contain additional 
information about the properties of LAEs and the inter- 
vening IGM. Due to the high scattering cross section of the 
Lya transition, modeling the diversity of Lya profiles ob- 
served in star forming galaxies requires radiative transfer 
calculations (Ahn etal. 2003; Hansen & Oh 2006; Ver- 
hamme, Schaerer, and Maaselli 2006; Tapken etal. 2007). 
Outflows are thought to decrease the intrinsic Lya atten- 
uation via frequency shifting, so their presence may make 
it easier to detect Lya emission. Since the Magellan LAEs 
have high SFRs, we interpret the line profiles in terms of 
the properties of galactic winds blown by the LAEs. 

We measure Lya lincwidths of 260 km s~^ and 
240 km s-i FWHM, respectively, for MSDM80-H3 and 
MSDM29.5-h5. The Lya fine in the Keck spectrum of 
MSDM71-5 is only 98 km s~^ wide, but the fine is broader 
in the low-resolution Magellan spectrum. We attribute the 
discrepancy to spatial variations across the object, which 
are weighted differently by seeing and slit alignment |f| We 
fitted very simple single-component and two-component 
models to the high-resolution line profiles. The one com- 
ponent model is a Gaussian with all flux blueward of the 
peak set to zero to mimic absorption by intergalactic gas. 
This kernel was convolved with another Gaussian that rep- 
resents the instrumental broadening. The 2-component 
model has an additional, weaker Gaussian added redward 
of the main line to simulate backscattering off an outflow. 

The lower right panels of Figures [31 [H and O show the 
resuhs. The proflle of MSDM29.5-I-5 is well fitted with 

^ The same object is detected by both observations because indeper 
spectra. 



a single component. The fitted instrumental smoothing, 
'--^ 110 km s^^ FWHM, is consistent with the resolution 
of the Keck-LRIS, follow-up spectrum. The model indi- 
cates a FWHM prior to attenuation of 400 km s^^ . The 
high-resolution spectrum of MSDM80-I-3 shows a red wing. 
For an instrumental resolution of 95 km s~^ , the fitted 
single-component model has a width of 300 km s~^ prior 
to attenuation. A spectrum with higher signal-to-noise 
ratio is needed to determine the reality of additional line 
components, such as a possible narrow component located 
300 km s""'^ redward of the main line. 

We fitted a model introduced by Dawson et al. (2002) 
and Hu et al. (2004) that is now in common use. Westra 
et al. (2005, 2006) fitted the model to two extremely lumi- 
nous LAEs at z — 5.7 from WFILAS. The two-component 
model for MSDM80-I-3 reminds one of their result for 
LAE Sll_5236, which has a second peak offset 400 km s^^ 
redward of Lya that is only 20-90 km s~^ wide. This 
model has also been fitted to the composite line profile 
of LAEs discovered serendipitously in the DEEP 2 redshift 
survey (Sawicki etal. 2007). The intrinsic width of the 
primary component — 350 km s^^ , and the component 
fitted to the red wing is offset ^ 420 km s~^ . These val- 
ues are similar to what we find for the Magellan LAEs 
even though the DEEP2 LAEs are more luminous (mean 
LLya = 6 X lO''^ erg s~^) and found at lower redshift (me- 
dian redshift 4.1). The substantial star formation rates 
suggest that winds are the likely source of the redshifted 
emission, but gas inflow can produce similar line profiles 
(Dijkstra etal. 2006). Strong winds would play a critical 
role facilitating the escape of ionizing photons (Fujita et al. 
2003). 

4.2. Number of LAEs Detected 

We confirm three LAEs from the four survey masks. 
The total unobscured slit area and the filter bandpass at 
FWHM sample a volume of 4.50 x 10'*/if(f Mpc^. The slit 
and the filter attenuate every LAE. The amount of at- 
tenuation is unknown for any individual source but is well 
understood statistically for the population. For any model 
of the luminosity function, we can accurately predict the 
attenuation of objects drawn from such a population and 
compute the number of sources detected given the com- 
pleteness limit of each of our four search observations. 

For parameterizations of the luminosity function {L^^,, 
a), where 

d<p{L) = (f>4L/L,)"e-^^^-'d{L/L^), (2) 

we computed the mean number of LAEs that our survey 
would detect. We calculated the volume from which a 
source is drawn using the slit width and redshift interval 
over which the source could have been detected. Because 
we start with the intrinsic luminosity (or flux), we can 
model the mean attenuation precisely. We discuss the re- 
sults obtained using our 90% completeness level. 

Figure (TU] illustrates the range of luminosity function 
parameters that yield ~ 3 detections. The shaded regions 
in Figure [TUl denote the allowed combinations of (0,,L*) 
by requiring that the model yield 1.37 to 5.92 LAEs (85% 
confidence limits, orange) or 31^218 LAEs (95% CL in yel- 

wavelength calibrations put the line at the same wavelength in both 
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low), where the hmits were computed from Poisson statis- 
tics (Gehrels 1986). The number density of LAEs in the 
Magellan MNS survey constrains the product well. 
The contours of equal number show how lower normaliza- 
tion values, require higher to yield a fixed num- 
ber of objects N. Since our data do not constrain the 
faint-end slope well, we show results for three assump- 
tions a — —1.2,-1.6, and — 2.0. In the limit a — —2.0, 
a normalization similar to that at z ^ 4.9 in the SDF, 
0* « 0.0055 Mpc-3 (Ouchi etal. 2003) requires L* lie 
in the range 1.78 x 10^^ to 5.97 x 10''^ ergs a'^ . As 
the faint-end slope becomes flatter, the characteristic lu- 
minosity must increase to yield 3 detected objects if the 
normalization is held fixed. 

Based on our confirmation rate to date, additional spec- 
troscopic follow-up will confirm, at most, a couple more 
LAEs in our emission-line survey. The number of remain- 
ing line sources for which a LAE identity has not been 
ruled out provides a strict upper limit. Our follow-up ob- 
servations were remarkably complete; only 4 objects in the 
lOH field, and 7 candidates in the 15H field, were not re- 
observed without the blocking filter. Some emission-line 
objects, however, were not-recovered in follow-up spec- 
troscopy (8 and 7 galaxies in the lOH and 15H fields, re- 
spectively). In the previous section, we argued that only 
two single-line objects in the lOH field were not likely fore- 
ground emitters. To compute an upper limit, however, 
we count all 9 single-line sources (including the confirmed 
LAE MSDM80+3). We found 7 single-line objects includ- 
ing 2 confirmed LAEs in the 15II field. The total num- 
ber of remaining LAE candidates is 42. For comparison, 
adopting the Shimasaku et al. luminosity functions with 
a = -2.0, -1.5, and - 1.0 yields 17, 14, and 12 LAEs, 
respectively, in our IMACS-MNS survey. That our upper 
limit is well above this confirms that most of our emission 
line objects that remain unidentified are not LAEs as we 
have argued. 

4.3. Cumulative LAE Luminosity Distribution 

Monte Carlo simulations of the luminosity distribution 
are the preferred method for determining the best luminos- 
ity function parameters. For a larger (or deeper) survey 
with more sources than our survey, we would estimate the 
luminosity function parameters by drawing sources ran- 
domly from a model of the luminosity function, assigning 
them random positions in the field and applying the detec- 
tion biases associated with the slit, filter, and data process- 
ing unique to MNS. The likelihood of the true parameters 
lying in any particular range could then be determined. 

With just three confirmed LAEs, however, we use a sim- 
ple variant of the 1/Vmax method (Felten 1976; Dawson 
et al. 2007) to estimate the cumulative luminosity distri- 
bution. We model the maximum volume, Vmax,i, within 
which each source could have been recovered. We use the 
slit width and redshift interval within which the object 
would have a flux greater than that of our survey's 90% 
completeness limit. Due to the high fraction of candi- 
dates followed up and recovered, the correction for unob- 
served candidates is small and omitted. Proceeding from 
the brightest source to the faintest, we add up the num- 
ber density of LAEs to obtain the cumulative luminosity 
distribution. Our best estimate is a lower bound in some 



sense because of the possible attenuation of each LAE. 
Table |4] shows this cumulative luminosity function. 

Figure [TT] compares this distribution (solid squares) to 
counts from other surveys at redshift 5.7. The open 
squares denote faint LAEs from Santos et al. (2004) that 
are gravitationally lensed by a foreground galaxy cluster. 
This ground-breaking survey of cluster caustics demon- 
strates the existence of less luminous LAEs; but the LAE 
number density is far more uncertain than the statistical 
error (shown). The survey volume is difficult to estimate 
and very small (i.e. subject to cosmic variance), so one 
may not wish to weight these points heavily when con- 
straining the luminosity function. The GRAPES produced 
three LAEs at z > 5; the density from Pirzkal et al. (2007), 
which includes lower redshift LAEs, is denoted by a cross 
in Figure [TT] The luminosities and redshifts of the LAEs 
discovered by Hu et al. (2004, open circles) with wide-field 
narrowband imaging provide the most direct comparison 
to our Magellan survey. The latter, like many of the imag- 
ing surveys, was carried out with the SuprimeCam mosaic 
CCD camera on Subaru but is unique among them; a high 
percentage of candidates have been confirmed spectroscop- 
ically. The number density of LAEs estimated from our 
Magellan MNS survey agrees very well with the result from 
the survey of Hu et al. , which sampled a volume compa- 
rable to Shimasaku et al. 

Number density estimates from the Large Area Lyman- 
Alpha survey still require significant extrapolation for un- 
confirmed LAEs at redshift 5.7 (Rhoads et al. 2001). Their 
best estimate (Malhotra & Rhoads 2004, open triangle) 
is, however, consistent with our new result. The data 
at redshift 5.7 will improve with the spectroscopic con- 
firmation of other large surveys with SuprimeCam (Mu- 
rayama etal. 2007; Ajiki etal. 2003, 2004, 2006). In 
Figure (TTJ we show LAEs, selected photometrically from 
SuprimeCam imaging of the Subaru Deep Field (SDF); 
Shimasaku et al. (2006) published spectroscopy for about 
one-third of the sample. For a faint-end powerlaw slope 
of a = —2.0, their fitted characteristic luminosity is 
L^ = l-Slog ^ 10'*'^ ^^S^ ^i*^ normahzation 0* = 
1.6lQ'y X 10^** Mpc^'^; and the integral of this luminosity 
function is plotted in Figure 1121 The LAE number den- 
sity of the Shimasaku et al. LAEs is 3-6 times higher than 
our result over the luminosity range logi = 'i'2.8tn.l- In- 
spection of the error bars in Figure [Tl] shows that only at 
the low end of this range are our results even marginally 
consistent. 

It is not obvious to us that cosmic variance is the ex- 
planation of this discrepancy. Stark, Loeb, & Ellis (2007) 
argue quantitatively that the variance in the Subaru Deep 
Field LAE Survey at redshift 5.7 is dominated by Poisson 
shot noise rather than cosmic variance. The area imaged, 
about 0.20 square degrees, is similar to that covered by 
both the LALA survey and Hu etal. (2004). The shal- 
low narrowband imaging survey, WFILAS (Westra et al. 
2006), covers a larger area. Only two of 7 LAE candidates 
have been confirmed in their 0.72 sq. degree field, but 
these two objects are very luminous; and we show their 
number density in Figure [TT] (hexagon) . To estimate the 
total variance in our survey, we used the on-line calculator 
described by Trenti & Stiavelli (2008). To simulate the ef- 
fect of the sparse sampling of the slits, we used a low value 



Spectroscopic Search for LAEs 



9 



for the survey completeness (10%) and the combined sur- 
vey areas in the lOH and 15H fields. For any halo filling 
fraction less than 70%, the total variance was dominated 
by the Poisson term rather than the clustering term. We 
conclude that cosmic variance is unlikely to dominate the 
uncertainty in our measurement of the LAE number den- 
sity. 

The Shimasaku ct al. result in the SDF is consistent with 
the WFILAS result considering their uncertainties, but the 
estimated values differ by a factor of five. Since several of 
these surveys probed large volumes, we see no reason to 
favor the newer Shimasaku et al. result until spectroscopic 
confirmation is completed. The combined data do not con- 
strain the break luminosity, L*, well in Figure [TT] 

4.4. The Lya Escape Fraction 

The luminosity distribution of i-dropouts has been mea- 
sured to very faint limits in the HUDF and a few other 
fields with deep HST/ACS observations (Bouwens etal. 
2007; Bouwens etal. 2006). Since we do not detect con- 
tinuum from our faintest LAEs, and spectra have not 
been acquired for all of the i-dropouts, it remains un- 
clear whether the Lya -selected sample is drawn from 
the i-dropout population. For purposes of illustration, 
we mapped the rest-frame UV luminosities, L^, of ob- 
jects drawn from the i-dropout luminosity function into 
line emission. We used the relation L{Lya ) = 1.47 x 
10^2 erg s"i/Lya(i,./102^erg s"! Hz"^), which assumes a 
timescale for star formation inherent to the mapping from 
UV luminosity to SFR given by Kennicutt (1998). Fig- 
ure [13 shows the cumulative Lya luminosity function of 
the i-dropouts. The three models are labeled by the Lya 
attenuation, fLya- For the no attenution example, this 
heuristic model overpredicts the number density of LAEs. 

The surprise is how easy it is to reconcile the num- 
ber density of LAEs and i-dropouts. Absorption by in- 
tergalactic gas significantly attenuates the Lya emission 
from z > 4 galaxies. For redshift ~ 6 galaxies, the SFR 
inferred from the UV luminosity is typically several times 
higher than that inferred from the Lya luminosity. The 
dashed curves in Figure [12] show that absorbing the blue 
half of the line profile, equivalent to setting fLya — 0.5, is 
one way to produce the LAE luminosity distribution from 
the i-dropout population. This particular model is not 
unique and may be a strong oversimplification of the rela- 
tion between these two galaxy populations. The ratio of 
UV to Lya luminosity presents enormous scatter among 
galaxies. The escape fraction f^ya could vary strongly 
with luminosity (or halo mass). 

Additionally, if the fraction of galaxies presenting Lya 
emission is less than unity, then the i-dropout luminos- 
ity function would need to be shifted in number besides 
star formation rate in order to model the LAE luminosity 
function. We know that only 25% of redshift 3 LBGs are 
LAEs (Shapley etal. 2003). The data for z ~ 6 galaxy 
samples appear to disagree. Shimasaku et al. (2006) and 
Kashikawa et al. (2006) argue for a rise in the fraction of 
galaxies with Lya emission to near unity at z 6, while 
Dow-Hygelund et al. (2007) found many i-dropout galax- 
ies without Lya emission. Modeling the effect of color cuts 
and equivalent width cuts on these samples will be required 

^ The recombination time depends on the clumping factor of the IGM. 



to clarify the Lya emission fraction. The correct physical 
mapping between the i-dropout and LAE populations is 
not currently known. Our comparison simply shows that 
a simple model with 50% Lya attenuation and a high frac- 
tion of LAEs among i-dropout galaxies is consistent with 
our estimated density of LAEs. That it severely overpre- 
dicts the counts from the lensed observations could mean 
a number of things: smaller Lya escape fractions from 
dwarfs, an overestimate of the number of faint i-dropouts, 
or an underestimate of the number density of lensed LAEs. 

4.5. Reionization Constraints 

The presence of Gunn-Peterson troughs in the spectra 
of z ~ 6 quasars (Becker et al. 2001; Djorgovski ct al. 2001) 
requires only a tiny neutral gas fraction in the IGM, about 
one in lO'' H atoms. The absence of strong evolution in the 
shape of the LAE luminosity function between redshift 5.7 
and 6.5 suggests little change in the neutral fraction lead- 
ing many to believe reionization was complete by redshift 
6.5 (Malhotra & Rhoads 2004; cf. Kashikawa etal. 2006). 
A protracted period of reionization from redshift 14 to 7 
would be consistent with the Thompson optical depth in- 
ferred from measurements of the cosmic microwave back- 
ground (Spergel etal. 2007). The time from redshift 7 
to 5.7 is 229 Myr, long enough that the intcrgalactic gas 
would recombine were it not bathed with Lyman contin- 
uum radiation(3 We therefore ask whether the LAE galaxy 
population discovered to date might supply these ionizing 
photons. 

To maintain IGM ionization at z = 5.7, the number of 
ionizing photons escaping from galaxies over the recombi- 
nation timescale must equal or exceed the number of hy- 
drogen atoms. The recombination rate is faster in denser 
regions of the IGM, so the required number of escaping 
photons increases with the dumpiness of the gas. The gas 
is smoothest on large scales and clumpier in the virialized 
regions where galaxies form. The gas in these halos ab- 
sorbs a high fraction of the ionizing photons (produced 
by galaxies), so the Lyman continuum escape fraction, 
fLyC, is lower on the larger scale of the IGM. For the 
regions of IGM penetrated by ionizing radiation, a clump- 
ing factor can be calculated that is the reciprocal of the 
volume filling factor of protons (or electrons) . The clump- 
ing factor is defined here in units of C =< rij^jj > /nj^jj 
(Madau etal. 1999); and at z '--^ 6 it is likely between 
1 and 10 (Furlanetto etal. 2006; Gnedin 2007). We pa- 
rameterize the clumping factor in units of C = 6 such 
that Cq = (1/6) < njjjj > /nHu- Following Madau 
etal. (1999), the critical ionization rate in the IGM is 
Nh = nH{0)/treciz) per unit comoving volume, which we 
can write as 

»„=5.2xlO"s-.Mp.-Ce(ii;y(2§)'. 

(3) 

The required production rate of H ionizing photons in 
galaxies is much higher because little Lyman continuum 
radiation escapes (Hurwitz et al. 1997; Leitherer et al. 
1995). Measurements at redshift 3 suggest absolute es- 
cape fractions fhyC^ 0.1 (Chen etal. 2007; Shapley etal. 
2006; Inoue etal. 2005; Fernandez-Soto etal. 2003; Gial- 

lowing Madau etal. (1995), Ir^^ = 290 Myr Cg"^ {tT^Y (i^)- 
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longo etal. 2002; Steidel etal. 2001). Numerical simula- 
tions suggest escape from the virialized halos of low mass 
galaxies is difficult (Gnedin et al. 2007) and that the re- 
sult is sensitive to the interplay between the blowout of a 
galactic wind and the star formation history (Fujita et al. 
2003). In the absence of any contribution from quasars, 
the minimum cosmic star formation rat^l required to keep 
the universe ionized at z=5.7 is 



p* 



0.02 Mq yr 

3 



LyCfl.l 



1 



z 



6.7 



-iMpc" 
0.047 



(4) 



where fhyC is the fraction of Lyman continuum photons 
escaping from galaxies, and fLyCfi.i means {fLyc/^-^)- 
The ratio C/ fLyc is not very sensitive to spatial scale, 
so the spatial dependencies of the clumping factor and es- 
cape fraction cancel each other out to some extent. 

To determine what fraction of this ionization require- 
ment the source population provides, we can add up all 
the LAEs, weighting by luminosity, in a given volume of 
space. Using Case B recombination conditions and Equa- 
tion[3l we find the critical Lya luminosity density required 
for full IGM ionization is 

3.0 X 10**° erg s^^Mpc"^ x 



-'Lya 



Ceil-O.lfLyCfl.i) 



Ihya 



0.5 



O ^2 \ 2 



.(5) 



jLyCfl.iJ V 6.7 J V 0.047, 
This equation highlights some obvious facts. If 
Lya radiation is highly attenuated (i.e. fhya « 1), then 
the observed Lya luminosity density can be low even 
though the observed galaxy population ionizes the IGM. 
Likewise, if Lyman continuum leakage from galaxies is very 
low (i.e. fhyC « l)i then even a high number density of 
galaxies will not completely ionize the IGM. To visualize 
the implications of the uncertain gas physics for the IGM 
ionization requirement, we collect all the factors in Equa- 
tion [5] describing the gas physics into a single parameter 

C = C(,{1 - 0.1fLyC,0.l)fLya,0.5fEyCfi.V (6) 

which is likely between 0.1 and 2. Although our simple 
modeling does not include possible variations in ^ among 
galaxies, we note that our survey limit is brighter than 
that expected at the theoretical mass cut-off for escape 
of Lyman continuum photons (Gnedin et al. 2007; Gnedin 
2007). 

We consider the luminosity density of the LAE popula- 
tion for two minimum luminosities: the Magellan MNS de- 
tection hmit of logLLyQ(erg s~^) — 42.57 and an extrapo- 
lation to logLLyaierg s~^) — 41.0. This lower value corre- 
sponds to a SFR « 0.1 Mq yr^^, a level of activity com- 
mon among dwarf starburst galaxies at low redshift. The 
line flux of such objects, F^ya ~ 3 x 10^^^ erg s^^ cm^^ 
at redshift 5.7, is well below current detection limits. 
For each luminosity function represented by the grid of 
Schechter function parameters in Figure 1101 we integrated 
the Lya luminosity density down to the two minimum 
luminosities. The total luminosity densities are shown in 
Figure [T31 Only those values within the (orange and yel- 
low) shaded region can be attained by a LAE population 
consistent with our number counts. 



The Lya luminosity density in Figure [T3] is shown rel- 
ative to that required for IGM ionization by Eqn. [5j Con- 
tours of equal luminosity density are labeled by the value of 
C required for ionization; blue and red contours correspond 
to a minimum luminosity rmnlog L Lya {erg/s) = 42.57 
and 41.0, respectively, for the lower integration limit. If 
spatial variations in Cq and fhyCOA do largely cancel 
out, then C, largely reflects the amount of Lya atten- 
uation, fhyafi.b- We consider values of Q from 0.1 to 
2.0 to be plausible. Comparison of the blue contours 
to the allowed (i.e. shaded) region shows that the de- 
tected population of LAEs can only ionize the IGM if 
logL*^iya(erg s~^) < 42.4. Following the arguments of 
§ 14. 4i the luminosity break for the i-dropouts, A/*_i6oo, 
maps to logL^xya — 42.60 for fLya = 0.5. Most previ- 
ous Lya studies have claimed L^^^ya greater than 42.4 
(Shimasaku etal. 2006; Malhotra & Rhoads 2004). 

For a population of LAEs with log L*,Lj/a (erg s-^)>42.5 
to ionize the IGM, galaxies much fainter than our detection 
limit must contribute to the ionizing radiation at z ~ 6. 
The blue contours do not intersect the allowed region in 
Figure [T^ for high values of \ogL^,Lya- Extending the lu- 
minosity function to very faint limits helps produce the 
stellar mass that appears to have been assembled by z ~ 5 
(Stark etal. 2007a). The flatter the faint-end slope of the 
LAE population, the harder it becomes to ionize the IGM 
when L*,Lj,Q is high. When logL^^^ya increases much be- 
yond 42.5 in Figure [151 the allowed (shaded) region is in- 
tersected only by the red contours with low values of C,. 
For large values of log L*, Lya (erg s~^) ~ 43, the parame- 
ter C, must be < 0.1 for complete ionization, much less if 
the faint-end slope is not steep. In our simple model, when 
\ogL^,Lya approaches 43, it becomes difficult for galaxies 
to ionize the IGM. 

A galaxy population described by lower values of L^,^Lya, 
however, easily ionizes the IGM. In Figure [T51 the inter- 
section of the blue contours with the shaded region indi- 
cates normalizations log(/)*,Lyc((Mpc~'^) greater than -2.2, 
-2.25, or -2.37 are required for faint-end slopes a = —2.0, - 
1.6, or -1.2, respectively. These high normalization models 
pull the luminosity break down to logL*xi;Q(erg s^^) = 
42.42, 42.38, and 42.35, respectively. These values of 
the break luminosity are below the detection limit of the 
Magellan MNS survey, so the models produce a larger 
Lya luminosity density per number of galaxies detected. 
If one accepts values of C as low as 0.1, then even the 
galaxy population detected by our Magellan MNS survey 
can keep the IGM ionized at z = 5.7. 

The volume of the Magellan MNS survey is not large 
enough to constrain well. One could as easily fit our 
points in Figure [11] with a power law instead of a Schechter 
function. To our surprise, however, we found that the 
other LAE surveys do not appear to constrain L.^, all that 
well. At redshift 5.7, a straight line can be drawn through 
the luminosity function data and error bars in Figure 1 of 
Malhotra & Rhoads (2004) , particularly when the number 
densities from lensing are excluded. The fitted distribution 
is better constrained in Figure 11 of Shimasaku et al. , but 
considerable leeway remains due to poor constraints at the 
bright end. We looked to lower redshift surveys for LAEs 



* Here, we have used an initial stellar mass function dN/dM oc M ^'^^ from 1.0 Mq to 100 M©. Continuous star formation at a rate of 
1.0 M0 yr~^ produces ionizing photons at a rate of Nh = 10^^"^ 8"^ (Kennicutt 1998; Leitherer etal. 1999). 
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to determine what values of might be reasonable. The 
MUSYC survey at 2; '--^ 3, in fact, favors significantly lower 
values, log L*,Lyaa (erg s^^) = 42.64^015 (Gronwall etal. 
2007); but inspection of their Figures 10 and 11 suggests 
the luminosity break is again not that well constrained 
due to the relatively shallow survey depth. It is possible 
that L^ Lya is lower at z = 5.7 than at redshift 3 since 
M^y may fade with lookback time from redshift 4 to 6 
(see Bouwens etal. 2007; cf. Iwata etal. 2007). 

We question whether these lower models have 

been ruled out by the data at this time. If not, then the 
population of LAEs detected could contribute a high frac- 
tion of the radiation needed to maintain the ionization of 
the IGM at redshift 5.7. If is large, then a large ex- 

trapolation of the luminosity function to fainter sources is 
required for LAEs to account for the ionizing photons. The 
faintest z = 5.7 galaxies recovered by the lensing surveys 
have unlensed luminosities of \ogLLya ^ 41.4, a factor of 
2.5 higher than the value of Lmm used to compute Fig- 
ure [131 Although there is tentative evidence for a steep 
faint-end slope (Stark etal. 2007b; Kashikawa etal. 2006), 
it seems likely that new technology will be required to de- 
termine definitively whether star-forming galaxies ionize 
the IGM at high redshift. 

5. EMISSION-LINE DETECTION OF HIGH-REDSHIFT 
GALAXIES 

We have demonstrated that multislit, narrowband spec- 
troscopy is a viable technique for finding LAEs at redshift 
5.7. Here we compare the efficiency of this sky-supression 
technique to other strategies using current instrumenta- 
tion. The sensitivity limits of next generation facilities to 
emission-lines are then discussed. Hydrogen Lya emis- 
sion will be the brightest line where the local gas kine- 
matics allow escape and the IGM neutral fraction allows 
significant transmission. Our discussion focuses on Lya 
detection over the period of reionization from 6 < 2: < 30 
but also applies to nearby lines like Hell 1640, CIV 1550, 
and NV 1240. 

The canonical figure-of-merit for galaxy surveys is the 
product AVl, where A is the collecting area of the telescope 
and Vl is the solid angle of sky observed at one time. The 
time required to survey a given area of sky to a particular 
magnitude is then inversely proportional to AU, for broad- 
band imaging. For emission-line surveys, the quantity to 
maximize, when sky emission is the dominant source of 
noise, is A^e/ {Y,\l^\9'^), where AA, 0^, and e represent 
spectral resolution, image size, and system throughput, 
respectively. For Lya observations from redshift 4 to 20 
(Ao = 0.6 — 2.6/zm), the surface brightness of the sky, 
is dominated by the zodiacal background in space. Its 
brightness is a strong function of distance from the eclip- 
tic plane and was measured with COBE (Hauser et al. 
1998; Kelsall etal. 1998). Atmospheric bands of molecular 
lines, mainly OH (hydroxyl and O2), dominate the broad- 
band background from the ground. The interline back- 
ground at high-dispersion is much lower than the broad- 
band mean, but its physical origin and intensity are not 
well understood. Maihara etal. (1993) suggest it is a few 
times brighter than the zodiacal background, but Bland- 
Hawthorne et al. (2004) argue that most of this emission is 
scattering and diffraction of the molecular lines. Here, we 



will assume the standard interline background from Mai- 
hara, and explore how the spectral resolution, image qual- 
ity, and system throughput inherent to various observing 
techniques affects both sensitivity and survey efficiency. 

Due to the high number of interlopers in many surveys 
for galaxies across the reionization era, this revised figure- 
of-merit does not completely describe the potential impact 
of a survey. The time required for redshift confirmation, 
not just that used to identify candidates must be consid- 
ered. In this regard, emission-line sensitivity is also of rel- 
evance for surveys employing continuum-break selection. 

5.1. Current Capabilities for Detecting Faint, 
Emission-Line Galaxies 

Emission-line spectroscopy with large, ground-based 
telescopes offers high sensitivity to objects with low SFRs. 
The HUDF observations with HST/ACS reached an ex- 
ceptional broadband depth. The luminosity function of i- 
dropouts reaches Muv,ab = —17.5 (Bouwens etal. 2007), 
a star formation rate of 0.6 Mq yr~^. The Lya emission 
from a subsample of these galaxies was detected in 30 hr 
integrations with Gemini (Stanway etal. 2006). Fainter 
LAEs clearly exist as shown by their detection with grav- 
itational lenses (EUis etal. 2001; Santos etal. 2004). 

Here we compare current capabilities of MNS, NB imag- 
ing, and tunable filters. Our arguments about the benefits 
of spectroscopy apply equally to observations with integral 
field units. The small fields of view of IFUs are inefficient 
for surveys at redshifts where the neutral fraction is low. 
At z ~ 10, however, the only detectable LAEs may be 
those within the ionized bubbles of the brightest sources. 
Furlanetto & Oh (2005) estimate physical sizes of 27 to 
227 kpc. Band- limited IFUs with diameters of 6-60" would 
be extremely effective for studying these galaxies. Slitless 
grism spectroscopy is never more sensitive than NB imag- 
ing because each pixel records background from the full 
transmitted bandpass. It has been used effectively from 
space owing to the lower background there and the ben- 
efits of low resolution spectra for identifying very faint 
galaxies (Xu etal. 2007; Pirzkal etal. 2007). 

5.1.1. Narrowband Imaging vs. Spectroscopy 

For any telescope and instrument for which emission- 
line detection is sky-noise limited, let us compare surveys 
carried out with MNS to those done with narrowband 
imaging. In a fixed amount of integration time t, spectro- 
scopic observations detect fainter emission lines because 
the bandpass of sky under the line, AA, is typically much 
less than the width of the filter, AAq. Maximum depth 
results from matching the spectroscopic resolution to the 
Lya linewidth. For a representative width of 150 km s~^ 
, the factor AA/AAq is approximately 4.lA/150A « 6~^; 
and the line flux reached with MNS is 

where quantities denoted by subscript zero represent nar- 
rowband imaging. Our IMACS pilot survey carried out 
with the facility 200-1 grism collected a 6.4 - 12. 8A band 
of sky under each emission line, so it did not reach max- 
imum depth. Many emission-line sources were unresolved 
by our observations, setting an upper limit on their size of 
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9< 0'.'75 FWHM. Better image quality would reduce sky 
noise further for both imaging and spectroscopy. The fac- 
tor of \/9{slit)/9t:, in Eqn.[7]a equals -^/I'/S/^o for our pilot 
survey. The gain appears larger for imaging because the 
sky aperture for MNS is fixed by the slit width in one spa- 
tial dimension, yielding a linear rather than a quadratic 
decreases in sky photons with seeing. For our pilot sur- 
vey, we cut slits twice as wide as the seeing to mitigate 
attenuation by the aperture. Hence, we think \/6/9q can 
be tuned to a factor y/2 in any seeing. The penalty in 
throughput, e, from the addition of the dispersing element 
is minor; for the 200-1 grism, the relative sensitivity is re- 
duced by just e/eo ~ 0.89^. We conclude that had we used 
our telescope time for NB imaging, our flux limit would 
be a factor of 3.0 to 4.3 times higher (i.e. less sensitive). 

When the total time required to search a given volume 
of space for emission-line sources is considered, the gain 
in sensitivity from MNS, relative to NB imaging, is off- 
set by the lower solid angle (per pointing) of the former. 
With the 200-1 grism, the clear aperture of the IMACS- 
MNS experiment subtends 55.3/656 = 8.4% of the active 
imaging field. In the time required to observe 12 masks, 
the integration time for direct narrowband imaging of a 
single field could be increased to to — 12t, increasing the 
line sensitivity to match that of MNS with IMACS. Nar- 
rowband imaging and MNS therefore reach comparable 
emission-line sensitivity when normalized by the volume 
searched per unit observing time. It is easy to show that 
this equivalence continues to hold as spectral resolution 
increases^ In a given amount of time, MNS reaches con- 
siderably fainter lines than NB imaging but selects from 
a smaller volume. When the objective is to optimize the 
yield of the faintest sources, the MNS technique is a better 
choice than NB imaging. 

The development of MNS will allow more sensitive sur- 
veys in the 8200A atmospheric window. Detector quan- 
tum efficiency, grism blaze, and spectral resolution can all 
be improved for IMACS observations, but installation of 
the new E2V CCDs will provide the single largest gain 
in throughput. In a deep ~ 25 hr survey, we expect to 
reach a line luminosity about three times fainter than the 
results presented here. The equivalent star formation rate 
at redshift 5.7 is just 0.79/£'J^ Mq yr'^ 

The biggest shortcoming of MNS is the time required 
to confirm candidates. Most of the emission-lines discov- 
ered with both MNS and NB imaging come from fore- 
ground galaxies. Narrowband imaging surveys parse the 
LAE candidate list using broadband colors and line equiv- 
alent width prior to spectroscopic confirmation. Broad- 
band surveys of the requisite area, however, do not go deep 
enough to identify the continuum counterparts of many 
emission-line galaxies discovered with MNS. Broadband 
depth will be less of a problem for differential imaging, 
another special technique under development for LAE sur- 
veys, because the field of view is contiguous and therefore 
of smaller extent for a given survey area. A good example 
is the 6'.8 by 6.'8 field of DAZzLE (Dark Ages redshift Lya 
Explorer, Horton et al. 2004) which subtends an area sim- 
ilar to that of our Venetian blinds mask and will be used 
to search for LAEs at at z '--^ 8 — 10. 



Possible solutions for rapid identification of LAEs from 
MNS data include resolved line profiles and equivalent 
widths. Recognizing the red skew of high-redshift Lya 
emission lines requires fairly high spectral resolution. The 
DEEP2 redshift survey provides a good example. Spectra 
were mostly taken with a 1200 line mm"-'^ grating. The ob- 
servations targetted z ^ 1 galaxies with multislit masks, 
but the "sky" regions of the slitlets subtend an area similar 
to our pilot survey. Sawicki et al. (2008) demonstrate that 
LAEs can be culled from the serendipitous detections of 
emission-line galaxies in the "sky" slits based on the shape 
of the line profile and the absence of additional emission 
lines. This second criterion is eliminated for MNS, and 
measurements of line asymmetry may become ambiguous 
at low signal-to-noise ratio. Deeper broadband imaging in 
even a single band longward of the emission line appears 
to be the fastest way to parse the emission-line candidates, 
as emission equivalent widths will normally be very high 
for Lya due to the combination of the {1 + z) boost and 
the potentially young, metal-poor stellar population. 

5.1.2. Tunable Filters 

Imaging observations with tunable filters can provide 
the same sensitivity and survey volume as MNS. The sky 
noise is comparable when the filter is cut down to the spec- 
tral resolution element of the MNS survey. The survey 
volumes will also be identical. While MNS slices the field 
spatially, it observes the entire redshift depth of the atmo- 
spheric window at once. Tunable filters image the whole 
field at once - one redshift slice at a time. The spatial and 
spectral slices are the same fraction of the volume for ob- 
servations with the same camera and spectral resolution. 
The advantage of imaging is exact positional information 
for the emission-line sources. The position - wavelength 
degeneracy inherent to MNS will present ambiguities when 
sources need to be identified in deep continuum images. 

We have taken data in the COSMOS field with the re- 
cently commissioned Maryland-Magellan Tunable Filter 
(MMTF) on IMACS. The highest resolution in the 8200A 
is about I2A, similar to our IMACS-MNS survey. We 
found a slightly higher throughput with the grism than the 
etalon, all other optical elements are the same. At redshift 
5.7 both search techniques are viable for future searches, 
with the MMTF offering better positional information and 
MNS offering factor of 1.5-2.0 better line sensitivity. Fur- 
ther to the red, the OH-free windows are narrower; and 
tunable filters offer a substantial advantage over MNS. 
Although the MMTF is not sensitive beyond the 9000 - 
9250A atmospheric window, other Fabry Perot spectro- 
graphs wiU be. The Flamingos-2 Tandem TF (F2T2) is 
an infrared-optimized scanning TF, which will be mounted 
inside the Flamingos-2 imaging spectrograph of Gemini 
south. The spectral resolution of AV ~ 375 km s~^ will 
deliver broader lines than optimal but will likely offer an 
efficient ground-based survey strategy at 2; > 7 in the near 
future. 

5.2. Relation to Next- Generation Facilities 

Spectroscopy with the next-generation of large ground- 
based telescopes will reach extraordinary sensitivity for 



^ The higher dispersion requires a larger separation between sUts. The increase in depth is a single pointing is cancelled by the increase in the 
number of pointings required to map the full aperture of an imaging survey. 
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emission-lines in the airglow band gaps. The gain from 
diffraction-hmited images will be as important as the in- 
creased aperture. Relative to our our IMACS MNS survey 
(subscript "P" for pilot), MNS observations with a next- 
generation ground based telescope of aperture D could 
reach 
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(8) 



where Cstrehi refers to the fraction of total light from an 
object in a circular aperature of diameter 9 and is related 
to the Strehl ratioH 

With adaptive optics and a 30 m aperture, line 
fluxes 21 times fainter than our IMACS observations 
become detectable. This sensitivity reaches SFRs of 
0.11 f£y^ M© yr~^ at redshift 5.7 in 10 hours. Such a pop- 
ulation may include the massive, young star clusters whose 
relics we see as the globular cluster population of the Milky 
Way today. 

Next generation surveys may actually detect fainter line 
fluxes from higher redshift objects. First, the image qual- 
ity delivered by adaptive optics will be better at longer 
wavelengths. Second, at high spectral resolution, the 
inter-hydroxyl-line background falls slowly with increas- 
ing wavelength between the 8200A bandpass and 1.9/x mF^ 
Ground-based line searches should become extremely sen- 
sitive to emission lines from redshift 5.7 to 15. Longward 
of 2/im , thermal emission causes the background to rise 
and observations with JWST will have a strong advantage. 

The ionized regions at redshift 7 to 10 will be mapped 
out by pointing these large telescopes at bright LAEs and 
finding fainter LAEs in the ionized bubbles. The Lya line 
will only be transmitted through the IGM when the source 
is surrounded by an ionized bubble (Furlanetto & Oh 
2005). An imaging tunable- filter, although not presently 
planned for first-light instrumentation on a 30m, could be 
the instrument of choice for such surveys. A more popular 
option for a facility instrument may be sets of deployable 
integral field units. A combined field-of-view exceeding 
6" X 6" could map the HII regions and provide spectral 
resolution or coverage for flagging interlopers. 

The extreme sensitivities attainable at the diffraction 
limit will only be realized if some Lya emission regions 
are small, < 50 pc. Some theoretical models do predict 
small sizes (Le Delliou etal. 2006), but Lya emission ha- 
los could also be quite extended. Empirical scaling re- 
lations for nearby HII regions, Rhii cr^ (Terlevich & 
Melnick 1981), suggest we should expect narrow lines of 
width < 40 km s"^ FWHM for small LAEs near the de- 
tection limit. The progenitors of the Milky Way's globular 
cluster systems could plausibly be detected as young star 
clusters. 

5.2.1. Ground-Based Emission- Line Sensitivity vs. 

JWST 



Observations from the James Webb Space Telescope will 
not be hindered by telluric- OH emission, which is pro- 
duced in a layer of atmosphere at roughly 87 km, or ab- 
sorption by molecules in earth's atmosphere, which pre- 
cludes observations over about 20% of the 1-5/im region. 
The background from space is dominated by sunlight scat- 
tered off zodiacal dust and is strongly dependent on eclip- 
tic latitude and longitude. Average values are about 10% 
as bright as the Maihara (1993) model for the interline 
background from the ground, although "deep fields" will 
presumably be carried out in regions of minimum zodi. 

The JWST instruments FGS-TFI and NIRSpec provide 
the highest sensitivity to line emission. The planned con- 
figurations of these instruments do not minimize sky back- 
ground. A single MEMS aperture on NIRSpec is much 
larger than a point source, and the spectral resolution of 
the Fabry Perot will not resolve line emission from galax- 
ies. These designs reflect the dominance of detector noise 
over shot noise from the sky for many observations. Be- 
cause JWST will be cooled to < 50 K, emission-line obser- 
vations with JWST will have their biggest advantage over 
ground-based line surveys at A > 2/xm . The background 
(in frequency units) should reach a minimum between 3 
and 4/xm due to the slow decline of the zodiacal back- 
ground with wavelength. In addition to extending studies 
of LAEs to redshift 30, the redder (and line-free) spec- 
tral coverage will be critical for studying hydrogen Balmer 
emission from objects toward the end of Reionization. 

Emission-line observations at z < 15 with JWST will 
not be restricted to OH-free bandpasses but are unlikely 
to be as sensitive as ground-based spectroscopy with a 
30 m class telescope. Using the background-limited case 
for scaling, the limiting flux from space relative to a 30 m 
telescope is 
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Assuming a relative sky brightness, Eao/SjvK ^ 10, 
ground-based spectroscopy could detect lines about a fac- 
tor of five fainter than JWST if the diffraction limit of a 
30 m telescope can be reached. The line sensitivity of these 
two telescopes is comparble at equivalent spectral resolu- 
tion and image quality; the extra collecting area roughly 
canceling the brighter sky. Even in the JWST era, large 
ground-based telescopes may lead the search to map the 
progression of reionization with LAEs. 

6. SUMMARY 

This paper presents the first LAEs discovered with nar- 
rowband, multislit spectroscopy. The technique succeeds 
on the 6.5 m Baade telescope because of the large detector 
on IMACS and its wide field-of-view. Our configuration 
reduced the sky noise by a factor of 10 relative to imaging 
through the same narrowband filter and provided spatial 
multiplexing 100 times that of a single longslit observation. 
The faintest confirmed LAE in the pilot survey has an ob- 

The Strehl ratio refers specifically to the ratio of peak intensity from a point source to the maximum intensity in a diffraction limited image 
with the same total flux. 

The Gemini ITC background model suggest the between-the-lines background in units of reaches a maximum in J band (around 
f.25/xm ) which is is about three times the 8200A background. Converting from frequency units to wavelength units cancels this out, 
Sa(1-25) _ „ .„ S^(1.25) 
E;^(0.82) " '-^E„(0.82)' 
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served line flux of F^ya = 1.4 ±0.2 x 10"^'' erg cm~^, 
which is near the estimated 90% completeness limit. 

The three redshift 5.7 LAEs have luminosities L^ya > 
5 X 10^^ erg s~^, which correspond to star formation rates 
greater than 5-7 Mq yr~^(for Salpeter IMF from 0.1 to 
100 M0 and solar metallicity). The fraction of intrinsic 
Lya emission absorbed by the galactic ISM likely varies 
with age and wind strength, while intervening IGM will 
absorb the Lya photons emitted blueward of the line 
center, contributing to the red skew of the line profile. 
The true star formation rates could easily be closer to 10- 
30 M0 yr~^. Winds with speeds of a few hundred km s^^ 
are common among local galaxies with star formation rates 
comparable to these LAEs and are consistent with the red 
wing of Lya line profile. The rest-frame ultraviolet contin- 
uum of these LAEs is fainter than the sensitivity limits of 
current wide-angle surveys. Future broadband measure- 
ments will constrain the age, mass, and dust content of 
these galaxies. Their relation to the i-dropout population 
is likely complex, although our best estimate of the cumu- 
lative Lya luminosity distribution is consistent with being 
drawn from the i-dropout population around redshift 6. 

Results from the IMACS MNS survey tightly constrain 
the product of the luminosity function parameters L*,lj,c( 
and 4'*,Lya- Nearly all candidates were followed up with 
both low- and high- resolution spectroscopy. The LAE 
population detected to date is unlikely to supply enough 
Lyman continuum photons to maintain the ionization 
of the IGM at redshift 5.7. For this detected popula- 
tion of LAEs to ionize the IGM, both a low value of 
the break luminosity ( and therefore a high normaliza- 
tion) and a combination of clumping, attenuation, and 
escape factors yielding C ~ 0.1, where ( = Cq{1 — 

0.1/Lj/C,o.i)/LyaQ,o.5/LyC.o.i' ^ould have to be realized. 
Determining whether galaxies ionize the IGM will clearly 
require deeper surveys that pin down the lumnosity func- 
tion parameters. 

Deeper emission-line surveys in the OH-free bandpass at 
8200A are possible with cither MNS or the MMTF at Mag- 
ellan. The former approach offers maximum depth while 



the tunable filter offers more accurate source positions. 
Both methods probe the same volume per pointing, which 
is much larger than the volume strongly lensed by a cluster. 
These techniques can be used on Magellan to reach lines 
2-4 times fainter than the results presented here. Such 
depth over a large volume should remove ambiguity about 
the value of I/»,Lj/a, particularly if complemented by shal- 
lower surveys with the volume to count LAEs brighter than 
\ogLLya{eTg s"i) = 43. 

Emission-line sensitivities will likely need to be im- 
proved further to find both the sources responsible for 
reionization and those that prevent the IGM from recom- 
bining at z '--^ 6. Finding these sources is a key science 
goal of the JWST mission. The next generation of ground- 
based telescopes with primary mirrors ~ 30 m in diameter 
(and adaptive optics) will offer even higher sensitivity to 
line emission shortward of 2//m , extending LAE searches 
to redshift 15 and detecting other lines like Hell 1640 at 
lower redshift. Our experiences with MNS indicate that 
isolated lines, lines with very large observed equivalent 
widths, and lines with red-skewed, asymmetric profiles are 
not exclusively Lya . Culling or exploiting the large num- 
ber of foreground emission-line galaxies from samples will 
have to be a key part of any experimental design. While 
identification of the Lyman break is one robust method, 
it obviates one perk of line searches, which is discovering 
objects below the detection limit of continuum surveys. 
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Foundation, the Alfred P. Sloan Foundation, and the Na- 
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Table 1 

Magellan Multislit, Narrowband Emission-Line Survey 



Date 


Mask 


RA 


DEC 


Slit PA 


Exposure 


Clouds/ 






(.J2()()()) 


(;j2()()()j 


(■^ ) 


(hrj 


Image Qualit \- 


2004 April 


COSMOS-0 


10:00:42.9 


+02:11:00 


90 


6.33 


Clear: 0.8-2.5" 


2004 April 


15H Field-0 


15:23:35.48 


-00:08:00.00 


135 


10.8 


Clear; 0.8-2.5" 


2005 March 


COSMOS-0.5 


10:00:42.9 


-1-02:11:00 


90 


10.0 


Clear; 0.6-0.8" 


2005 May 


15H Field-0.5 


15:23:35.48 


-00:08:00.00 


135 


6.75 


Thin Clouds; 0.6-0.8" 
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Table 2 

Candidate Yield and Spectroscopic Classification 



(1) 




(2) 


(3) 


(4) 


(5) 


(6) 


Survey 




Confirmed 


Confirmed 


Foreground 


No Follow 


Not 






LAEs 


Single Line 


Objects 


Up 


Recovered 


lOH 2004 - 


Line Only 


1 


7 


20 


3 


5 


lOH 2005 - 


Line Only 





2 


27 


1 


3 


15H 2004 - 


Line Only 


1 


6 


25 


6 


4 


15H 2005 - 


Line Only 


1 


1 


19 


1 


3 


lOH 2004 - 


Line + Cent. 








18 


16 


1 


lOH 2005 - 


Line + Cent. 








20 


12 


4 


15H 2004 - 


Line + Cont. 








21 


11 


8 


15H 2005 - 


Line + Cont. 








19 





1 



Note. — Cols - (1) Survey name and description of the continuum flux in the NB8190 band. (2) Number of confirmed LAEs. (3) Number 
of LAE candidates confirmed to be single-line sources. (4) Number of lines in 8190 band identifed as foreground objects. (5) Number of 
emission-line objects not followed up. (6) The Number of emission-line objects not redetected. Note that exposure times of follow vary. 
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Table 3 
LAE Candidates 



(1) 


(2) 






(3) 


(4) 




(5) 


(6) 


Name 


RA 






DEC 


F(Lya) 

( erg s ^ cm 




Wo 


Ao 




(J2000.0) 


(J2000.0) 




(A) 


(A) 


MSDM80+3 


10:00: 


:30, 


.479 


02:17:14.80 


14± 2 X 10" 


18 


> 244 


8141 


MSDM64-5 


10:01: 


:12, 


.769 


02:13:14.80 


63±8 X 10- 


18 


> 1082 


8141 


MSDM13-3 


10:00: 


:54, 


.458 


02:00:30.67 


14± 2 X 10" 


18 


> 240 


8169 


MSDM50-7 


10:01: 


:23 


.675 


02:09:45.15 


16±4 X 10" 


18 


> 281 


8122 


MSDM52+8 


09:59: 


:53, 


.266 


02:10:15.11 


25 X 10-1^ 




421 


8179 


MSDM30.5-8 


10:01: 


:29, 


.407 


02:04:53.60 


26±8 X 10- 


18 


> 447 


8173 


MSDM17.5+8 


10:00: 


:08, 


.176 


02:01:38.46 


43 ± 9 X 10" 


18 


> 739 


8129 


MSDM 17-1 


10:00: 


:45, 


.109 


02:01:30.41 


32 ± 2 X 10" 


18 


> 555 


8162 


MSDM29.5+5 


15:22: 


:57, 


.880 


-00:07:36.34 


18±2 X 10" 


18 




8148 


MSDM71-5 


15:24: 


:08, 


.929 


-00:10:43.06 


22±6 X 10- 


18 


> 61 


8130 


MSDM89-2 


15:24: 


:03, 


.935 


-00:03:06.53 


14±5 X 10- 


18 


> 39 


8200 


MSDM57+4 


15:23: 


:19, 


.653 


-00:03:20.27 


21 ±6 X 10" 


18 


> 58 


8178 


MSDM36+5 


15:23: 


:03, 


.246 


-00:06:38.71 


17±8 X lo- 


18 


> 47 


8135 


MSDM54.5-1 


15:23: 


:37, 


.234 


-00:08:36.91 


ll ±6 X 10- 


18 


> 31 


8222 


MSDM55-3 


15:23: 


:48, 


.230 


-00:11:11.24 


4.9± 5 X 10" 


-18 


> 14 


8170 


MSDM94-4 


15:24: 


:15, 


.670 


-00:04:16.34 


9.5± 7 X 10- 


-18 


> 26 


8144 


MSDM70+1 


15:23: 


:44, 


.420 


-00:04:56.60 


9.4±5 X 10" 


-18 


> 26 


8227 



NOTES - In addition to confirmed LAEs MSDM80-I-3, MSDM29.5-F5, and MSDM71-5, we list candidates 
presenting a single line in our low resolution follow-up spectra. (col 1) - Object Name. Format is 
MSDM[Slit Number]. [Field Identifier] [±Slit Segment] . (col 2, 3) - Coordinates. See § [O for full description of un- 
certainties, (col 4) - Line flux measured in search image, (col 5) - Measured emission equivalent width in the observed 
frame. An upper limit on the continuum flux density of 5.8 x 10"^° and 3.6 x 10"^^ erg s~^ cm~^ is assumed for the 
lOH Subaru z imaging and the 15H LCIRS z imaging, respectively. An entry of nodata indicates that this imaging does 
not cover the position of the object, (col 6) - Wavelength measured in follow-up spectrum. The uncertainty is dominated 
by the placement of the object in the 1'.'5 slit. One half the slit width corresponds to a wavelength uncertainty of about 
7'.' 5. 



Spectroscopic Search for LAEs 



19 



Table 4 

Cumulative Luminosity Distribution of LAEs 



F 




L 


V 


0(> L) 95%CL 


(^(> L) 84.1% CL 


( erg s"^ 


cm-^) 


(erg s-i) 


(/ifo' Mpc-3) 


mo Mpc-3) 


(/ifo Mpc-3) 


2.2 X 10" 


-17 


7.8 X 10*2 


4.429 X 10* 


2.26l;^i' X 10"^ 


2.26i5'-9 X 10-5 


1.8 X 10" 


-17 


6.4 X 


3.904 X 10'' 


4.79+^'\ X 10-5 


4.79t^'^ X 10-5 


1.4 X 10" 


-17 


5.0 X lO''^ 


2.402 X W 


8.95t^"i X 10-5 


8.95t4^4 X 10-5 



Notes: Col(l) Estimated line flux of LAE. Col(2) - Estimated Lya luminosity. Col(3) Survey volume within which the 
object could have been detected. Col(4,5) Cumulative number density of LAEs computed from (f>(L > Li) = TijUj ~ 
T,j{l/Vm,j), where the sum is over all galaxies j with Lj > Li. We divide the Poisson error on 1 object by the volume 
Vj to get the error on each nj = 1/Vj term. Following Gehrels (1986), the 84.1% (95%) confidence range for 1 detection 
is 0.173 to 3.300 (0.0513 to 4.744) objects, respectively. And the total error is Sep = S{N = l)[Sj(l/Kij)^]"-5, where 
(5(iV = 1) takes the values +3.744 and -0.9487 objects (in column 4) and +2.300 and -0.827 objects (in column 5). 
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Fig. 2. — IMACS slit mask used for the spectroscopic survey. Coordinates denote positions in the focal plane. The I'.'S wide slits subtend 
a geometric area of 55.3 □ per pointing within the 27.'50 diameter field of the camera. 
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Fig. 3. — LAE MSDM29.5+5. The zeropoint of the velocity scale is arbitrary, defined with zero at maximum intensity for convenience. Tics 
mark the separations of lines that commonly appear in our spectra of foreground interlopers, (a): IMACS 200-1 f/2 spectrum taken in 2005 
June. The line is narrow compared to the profile of a sky line, shown by the light-blue, dotted line, (b): Spectrum from (a) on a wavelength 
scale. Also shown is an IMACS 150-1 f/2 spectrum taken taken in 2006 June (thick, blue line), (c): Higher resolution, Keck/LRIS spectrum. 
The line is asymmetric and skewed to the red. (d): Keck/LRIS spectrum and fit. The model is a Gaussian (400 km s-^ FWHM) with flux 
blueward of the peak set to zero and then convolved with another Gaussian representing instrumental broadening (110 km ). 
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MSDM 




V (km/s) 



Fig. 4. — LAE MSDM80+3. Velocity scale, tic marks, and arc line as in Figure m a,b: IMACS f/2 spectra from 2005 May (thin, white 
line) and 2007 February (thick, blue line), c; IMACS f/4 spectrum, (d): IMACS f/4 spectrum and fit. The model is a Gaussian (330 km s"! 
FWHM) with flux blueward of the peak set to zero and a 150 km s~^ wide Gaussian 330 km redward of the first. These components 
were convolved with another Gaussian representing instrumental broadening (150 km s~^ ). 
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III 

- MSDM71-5 ® 




1 1 


1 

ki +[NII] 


III] 


- 













1 



-6000 -4000 -3000 2000 4000 6000 8130 8140 

V (km/s) 




V (km/s) 



Fig. 5. — LAE MSDM71-5. Velocity scale, tic marks, and arc line as in Figure|3] Top: IMACS f/2 spectra with the 200-1 grism in 2004 July 
(thin, white line) and the 150-1 grism in 2006 June (thick, blue line). Bottom: IMACS f/4 spectrum. The modelis a Gaussian (150 km s~^ 
FWHM) with flux blueward of the peak set to zero and then smoothed by an instrumental response of width 75 km . 
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Fig. 6. — Low-resolution (f/2) spectra of lOH Field LAE candidates. The objects shown have a single line in the observed-frame 5000-9000A 
bandpass. This inset shows no lines are seen at the spacing (tic marks) of the foreground contaminants [OIII]5007,4959 or Ha -|- [NIIJ6584. 
Comparison to the profile (dotted line) of a source that fills the slit shows the lines are unresolved, and higher resolution spectroscopy is 
required to show line asymmetry and/or identify the [Oil] 3726,29 doublet (small tics). 
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Fic. 7. Low-resolution IMACS spectra of 15H Field LAE candidates. The objects shown have a single line in the observed-frame 
5000-9000A bandpass. This inset shows no lines are seen at the spacing (tic marks) of the foreground contaminants [OIII] 5007,4959 or 
Ha + [Nil] 6584. Comparison to the profile (dotted line) of a source that fills the slit shows the lines are unresolved, and higher resolution 
spectroscopy is required to show line asymmetry and/or identify the [OH] 3726,29 doublet (small tics). 
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Line Flux (ergs/s/cm ) 



Fig. 8. — Completeness of emission-line surveys for two limiting sources sizes. For the 2005 search in the lOH Field, the line flux corresponding 
to 50% completeness is 6.6 - 7.2 X 10~l* erg s"! cm'^ for sources between 0''45 and 0'.'90 FWHM. 
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1.5 arcsec slit, posn from centre 




-2-10 1 2 

position (arcsec) 

Fig. 9. — Attenuation of flux vs. offset from center of 1'.'5 wide slit. Slit losses are shown for seeing-convovled sources with Gaussian 
brightness profiles of FWHM 0.5 to 1'.'2 in steps from O'.'l. The red curve represents the worst effective image quality in our survey. 
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42 42.5 43 




Fig. 10. — Mean number of LAEs detected in our experiment for different values of the intrinsic luminosity function parameters. From top 
to bottom, the panels show faint-end slope a = —1.2, 1 — .6, and — 2.0. The y-axis is the normalization of the luminosity function, 0*, The 
X-axis is the characteristic luminosity _L». The orange and yellow shading marks the 84.1% and 95% confidence regions for dectecting 3 LAEs. 
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95% CL 

84% CL (solid line) 



cm. 0001 
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■ Martin el al. 3008 
XPirzkal el aL 2007 (z - 
OWestra et al. 3006 
(^Shimasaku et al. 2006 
OHu et al. 3004 
□ Santos et al. 2004 
ARhoads et al. 3001 



4.9) 




{) 



log L (ergs s ^) 



Fig. 11. — Cumulative number density of LAEs. Our best estimate is shown by the solid squares; the errorbars are discussed in the caption 
to Table |4] We compare number densities from the following narrowband imaging surveys: LALA 2001 (triangle), Hu etal. 2004 (circle), 
Westra etal. 2006 (hexagon), and Shimasaku etal. 2006 (pentagons). Spectroscopic results are shown from GRAPES (Pirzkal etal. 2007, 
crosse) and a program that positioned longslits along cluster caustics (Santos etal. 2004, open squares). 
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Fig. 12. — Comparison of cumulative number density of z = 5.7 LAEs and i-dropouts. Blue curves show the fitted luminosity distribution of 
LAEs from Shimasaku et al. (2006) with faint-end slope a = —2.0, —l.5{solid), and —1.0. Black curves show the adopted fit to the distribution 
of i-dropouts from Bouwens etal. (2007); for purposes of illustration, we adopt constant hya escape fractions of /jy^ = 1.0,0.5, and 0.3. 
The dashed, black curve illustrates the magnitude of some systematic uncertainties in the luminosity distribution fit; and solid circles show 
the i-dropout counts, converted to Lya luminosity using /jy^, = 1.0. 
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42 
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Fig. 13. — Total hya luminosity density for the 3 values of the faint-end slope. Two lower limits are illustrated: (1) our survey detection 
limit of logL = 42.57 erg s~^ (blue contours), and (2) an extrapolation to logL = 41.00 erg s~^ (red contours), which corresponds to dwarf 
galaxies with SFR = 0-095 f^^^. The critical number of sources required to keep the IGM ionized is marked in units of the parameters 

Cs{l — 0.1fLyC,o.i)fLyaa,o.5/fLyC,0.l = l-O- With the lower limit in (1), it appears unlikely that enough galaxies have been discovered to 
keep the IGM ionized at z=5.7. It appears, however, from (2) that an extension of this population to lower luminosities, even with a shallow 
faint-end slope like a = —1.2, could easily produce the missing ionizing photons. 



